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Abstract 

We analyse new exact enumeration data for self-avoiding polygons, counted by 
perimeter and area on the square, triangular and hexagonal lattices. In extending 
earlier analyses, we focus on the perimeter moments in the vicinity of the bicriti- 
cal point. We also consider the shape of the critical curve near the bicritical point, 
which describes the crossover to the branched polymer phase. Our recently conjec- 
tured expression for the scaling function of rooted self-avoiding polygons is further 
supported. For (unrooted) self- avoiding polygons, the analysis reveals the presence 
of an additional additive term with a new universal amplitude. We conjecture the 
exact value of this amplitude. 

1 Introduction 

The model of (planar) self-avoiding polygons fTj ^] is an important unsolved model of 
statistical physics. Some progress has been made in recent years. Results from exactly 
solvable polygon models led to a prediction of the scaling function of self-avoiding polygons, 
counted by perimeter and area [22 E]- On the other hand, the theory of stochastic processes 
provided new insight into the problem by relating it to the so-called Schramm-Loewner 
evolution ^H]- In this article, we will further test the predictions implied by the scaling 
function conjecture J221I21 and extend it. 

Let Pm,n denote the number of self-avoiding polygons (SAP) of perimeter m and area 
n on a given lattice. In this article, we will consider self-avoiding polygons on the square, 
hexagonal and triangular lattices. Denote the perimeter and area generating function by 



G{x, q) = YlimnPm,nX'^'f- The function G{x, 1) is called the perimeter generating function 
of the model, its radius of convergence is denoted by Xc- For SAP, we have for the singular 
part of the perimeter generating function G'-'^™^-*(a;, 1) ~ A{l—x/xcY~^ as a; — > Xc, with a = 
1/2 being universally accepted, though not rigorously proved. (All limits appearing in this 
paper will be taken from below.) Until section 5, we will concentrate on rooted self-avoiding 
polygons, whose perimeter and area generating function is G^'^\x,q) = 'Yl,mnPrn\nx"^q^, 
where pm,n = fnpm^n- We thus have G^'^\x,q) = x-^G{x,q) and G'^'')''^'**"3)(x, 1) ~ B{1 — 

ryt / rf 1 -'- *-*■ Q Q /y \ If 

The phase diagram of SAP, enumerated by perimeter and area, appears to have first 
been discussed in [9]. There is a phase boundary in the region q < 1, x > Xc, terminating 
in a bicritical point at (xc, 1). If g = 1, typical polygons are extended, whereas for g < 1, 
typical polygons try to minimise their area. This particular type of phase transition is 
also called a collapse transition, and the phase g < 1 is called the branched polymer 
phase. Indeed, polygons of minimal area may be viewed as branched polymers. The phase 
boundary is characterised by a logarithmic singularity, when approached from below. This 
particular feature was first found by studying the area generating function of SAP, G{l,q), 
which was found (numerically) to have a singularity of the form G{1, q) ~ y41og(l — q) as 
g — > 1 [?]. The point (x, g) = (xc, 1) is a bicritical point where a generic scaling form of 
the perimeter and area generating function is expected to hold jU] . The singular behaviour 
about a bicritical point is generally expected to be of the form 

G(x, g) ~ G(-^)(x, g) + (1 - qfH ( (f^) (^, <i) - (^c 1), (1.1) 

where G^^'^^\x, q) denotes the regular part of G{x, q) at (x, g) = (xc, 1), and H{s) is called 
the scaling function with critical exponents and 0, and s = Tff-^- We stress that 
there are counter-examples known where such a scaling form is not valid, for example, in 
the simple model of rectangles [12]. For staircase polygons jTHI, the behaviour ()1.1|) has 
been proved. The exponents are 6 = 1/3, = 2/3, and the scaling function of staircase 
polygons is the logarithmic derivative of an Airy function. The phase diagram of staircase 
polygons is similar to that of SAP. There is also a phase boundary in the region g < 1, 
X > Xc, terminating in a bicritical point at (xc, 1). The phase boundary in that case is 
characterised by a simple pole, when approached from below, while at the bicritical point 
we have a branch-point singularity (j^**"'^'*(xc, g) ~ -8(1 — g)^''^ as g — >■ 1, as follows from 
^^. Interestingly, rooted SAPs display the same singularity structure as the staircase 
polygon model. This led to the question whether the scaling functions might be the same 



Any conjectured form for the singular behaviour of G*^''^(x, g) can most appropriately 
be tested by comparing predicted moments to those calculated numerically. In [1^1 122] 
we made such a comparison with the area moments, which led us to conjecture the exact 
form of the scaling function, thereby answering the previous question in the affirmative. 
We review this calculation below. The validity of the scaling function conjecture leads, 
in addition, to predictions of the leading singular behaviour of the perimeter moments, as 
explained in section 3. Checking this behaviour thus yields a further test of the scaling 



assumption. This was done in P21 only for the moment of order zero, i.e., the (bicritical) 
area generating function (?'•'"•' (xc, q) as q ^ 1. Here we present, for the first time, a detailed 
numerical analysis for the higher moments, which turns out to be numerically much more 
difficult than the analysis of the area moments. 

After establishing agreement with the predictions of the values of the first 10 perimeter 
moments in section 4, we consider the scaling functions for (unrooted) SAP, obtained by 
integration of the rooted SAP scaling function. The "constant" of integration, which must 
be a function solely of q, in order that its derivative with respect to x vanishes, turns out to 
dominate the behaviour of the scaling function as g — > 1. We argue for a particular form of 
this term, and then numerical testing reveals an unexpected amplitude universality across 
the three lattices we study. Based on our experience with other exact amplitudes, we 
conjecture the exact value of this universal amplitude. More precisely, we show below that 
at the bicritical point the behaviour is G^^^'^^\xc, q) ~ A{1 — q) log(l — g) as g — ;> 1, where 
we conjecture the exact value of the amplitude A. Our findings imply that the scaling form 
(II. ip cannot hold for (unrooted) SAP. We suggest a modified form below, see (15. 4|) . (|5.6p 
and dSHj). 

In the last section, we analyse the shape of the critical curve near the bicritical point, 
which describes the crossover to the branched polymer phase. We find the prediction 
from the scaling function conjecture satisfied, within numerical accuracy. We also analyse 
the behaviour of the critical curve as g — > 0. Numerical techniques are explained in an 
appendix. 

2 Area moments for rooted SAP 

The factorial area moment generating functions g)^ {x) of rooted self-avoiding polygons are 
defined by 

(-1)' 



«''<^> = Ta?'^"<^-') 



q=l «'• 



Y.^n),p%x^, (2.1) 



where (a)^ = a{a — 1) ■ ■ ■ {a — k + 1). Previous numerical analyses [1^1 123 based on exact 
enumeration data provide strong evidence for the asymptotic form 

with exponents 7^ = (3 A; — l)/2. We incorporate the coefficients and exponents into the 
function F^^\s) defined by 



F'''^') = Y.^.- (2-3) 



fe=0 



Slh 



At this stage, the function F^^\s) should be viewed as some generating function for the 
numbers fk, with s being an undetermined variable. We will argue below that, given the 
validity of the scaling assumption (jl.lj) . the function F^^\s) is the scaling function H^'^\s) 



Amplitude Square Hexagonal Triangular 

Eq 0.56230130(2) 1.27192995(10) 0.2639393(2) 

El 0.0795773(2) 0.0795779(5) 0.0795765(10) 

/o -1.61880474(6) -3.0645083(3) -1.906228(1) 

/i -0.01508198(4) -0.0215332(3) -0.0191714(2) 

Table 1: Values of area moment amplitudes for SAPs taken from 

of the model. ^ (In fact, F^'^\s) is then the generating function of the perimeter moment 
amplitudes e^, which are defined in ()3.2|) ). 

In |22], we tested the conjecture that the function F^'''\s) of rooted self-avoiding poly- 
gons is given by 

^.)W = _4/.AlogAi((Ay'%^ ,2.4) 

by comparing numerical estimates of fk with estimates that follow from ()2.4|) . The con- 
stants /o and /i have been numerically determined previously to great accuracy |21]. 
We have /o = — ^"^ and /i = — ^^^, where Xc is the radius of convergence of the 
perimeter generating function, given by Xc = 0.379052277757(5) on the square lattice. 



Xc = 0.2409175745(3) on the triangular lattice, and Xc = l/v 2 -|- \/2 on the hexagonal 
lattice PHI- The constant a is defined such that Pm,n is nonzero if m is divisible by a. 
Thus (7 = 2 for the square and hexagonal lattices and a = 1 for the triangular lattice. 
Estimates of the amplitudes Eq and Ei, taken from j21j, are given in Tabled The value 
El = l/(47r) = 0.07957747. . . has been derived in jl], using field theoretic arguments. The 
conjecture fl2.4|) follows from the assumption that rooted self-avoiding polygons behave 
asymptotically like models whose perimeter and area generating function are described by 
a g-algebraic functional equation of arbitrary degree with a square-root singularity as the 
dominant singularity of the perimeter generating function jH 122 I2B] • This class includes a 
number of exactly solved polygons models such as staircase polygons, column-convex poly- 
gons, and bar-graph polygons [201 112] • For models within this class, the coefficients fk and 
exponents •jk from ()2.2|) have been exphcitly calculated, leading to the expression ()2.4|) . It 
is interesting to note that the distribution of area in the limit of large perimeter, which can 
be extracted from ()2.4j) . see also '25', Appendix], is given by the Airy distribution, which 
appears in a number of related contexts [lti|2.5|. 

Assuming the validity of the scaling form (jl.lj) . the singular behaviour of the area mo- 
ment generating functions ()2.2j) determines the critical exponents and the scaling function 
H^^\s). Taking the limit g ^ 1 in (jl.ll) . we infer from ()2.2j) . leaving 7^ unspecified for the 

moment, that 

k-9 , , 

7fc = ^. (2.5) 

We also find that the asymptotic expansion of the scaling function H^^\s) is given by 
^The superscript ^'"^ indicates the relation to rooted SAP. 



H^'^\s) = F^'^\s), where F^^\s) is defined in ()2.3p . Thus, for rooted SAP, the assumption 
of the scahng form (ll.lj) . together with the result ()2.4|) for the area moments, leads to 
exponents 9 = 1/3 and (p = 2/3 and to the scaling function H^^'{s) = F^^'{s), where 
F(^)(s) is given in (jT^ . 

3 Perimeter raoments 

The factorial perimeter moment generating functions h}^ (g) of rooted self-avoiding poly- 
gons at the bicritical point are defined by 

ht\Q) = ^ ^^^'^(^'^) = ^E("^)^PSn^^^", (3-1) 

^—^c m,n 

The scaling assumption ()l.ip leads to a prediction for the behaviour of the perimeter 
and area generating function in the limit x — >■ Xc- The leading singular behaviour of the 
perimeter moments is given by 

hV'^'''''^^) - J^^ (^-1)> (3-2) 

with exponents j3k = kef) — 6, where the singular amplitudes e^ appear in the expansion of 
the scaling function H^'''\s) = F^'^'{s) about the origin. 



oo 



F(^){s) = Y,eks'. (3.3) 

fc=0 

We can readily derive formulae for the expansion coefficients e^ in ()3.3p . Note that the 
scaling function (j2.4p satisfies the Riccati equation 

fW(.)2-4AFW'(s)-/2 = 0. (3.4) 

Inserting the form ()3.3|) into ()3.4|) . we obtain for the numbers e„ the expression 

/ . N (2n+2)/3 

en = bnfiij-) ' (3.5) 

where the constants 6„ are defined by the quadratic recursion 

n-l 

Anbn + Sn-1,1 = ^ bkbn-i-k (n > 0). (3.6) 

fc=0 

The constants 6„ are polynomials of degree n + 1 in bg. The value of bg, remaining unde- 
termined by (j3.4j) . can be extracted from the limit s — i> in (j2.4j) as 

6„ = -^— (3.7) 

5 



It is interesting to note that the perimeter moments are related to the area moments 
of negative order |Tni Eqn. (37)]. It follows from ()3.5p that the amplitude combinations 
ekei''eQ~^ = hkhi^hQ~^ are independent of the amplitudes /o and /i. The first few combi- 



nations are 



v^TT-^ 



6261% = 1-—— 2-^, 6361^6^ = 1 

_4 3 _ ^ 10 V^TT^ -5 4_T 

6461 Cq — i 8irf2\6' ^5^1 ^Q — ^ 



8 v^7r3 






8ir(|)6' 






4 VSvr^ 
27r(|)6 ^ 


16 vr^ 
1215r(|)i2 



(3.8) 



In the next section we compare these predictions with the numerical values obtained 
from our new enumeration data. 

4 Perimeter raoment analysis 

We have generated data for self-avoiding polygons, counted by perimeter and area on the 
square, hexagonal and triangular lattices, using the finite lattice method. In particular, 
we determined the numbers Prn,n for n < 50 (square lattice), n < 40 (hexagonal lattice) 
and n < 60 (triangular lattice), for all relevant perimeter lengths. The algorithms used 
in our SAP enumerations are based on the finite-lattice method devised by Enting [^ in 
his pioneering work on the enumeration of polygons on the square lattice. Details of the 
algorithms used to enumerate SAPs on the hexagonal and triangular lattices can be found 
in jH] and |HI, respectively. A major enhancement, resulting in exponentially more efficient 
algorithms, is described in some detail in |^ while recent work on parallel versions can 
be found in ^Sj. All of the algorithms described in these papers are for enumerations by 
perimeter, but the generalisation to include area is straightforward. The calculations were 
performed on the server cluster of the Australian Partnership for Advanced Computing 
(APAC). The calculations for the square lattice required up to 14Gb of memory, and were 
performed on up to 16 processors using a total of just under 2000 CPU hours. Comparable 
computational resources were required for the hexagonal and triangular lattices. 

We first checked the prediction for the exponents (5k = 2fc/3 — 1/3 defined below ()3.2p . 
using first order differential approximants |llj. Then, we estimated the amplitudes eo and 
ei by a direct fit of the data to the expected asymptotic form, as explained in the appendix. 
Using the notation of the appendix, we fitted with exponents of the form ai = {i + l)/3. 
In the data analysis, we had 1 < i < Mq, where 2 < Mq < 4. 

The particular choice of exponents ttj arises from the numerically well established be- 
haviour of the area moment generating function 

1=0 ^ '^ ' 



with exponents 



^k,i = {3k-l- l)/2, (4.2) 



where fk,i = and f 0,21+1 = 0, see |22I- If, in generalising ()1.H) . a scaling behaviour of the 
form 



1=0 ^ 



qY 



'"' "" ' Xg)^(a;e,l). (4.3) 



is assumed with exponents 61+1 > 61, then by the arguments of the last section, the limit 
q —^ 1 constrains the exponents to 

lk,i = -^. (4.4) 

Comparison with ()4.2|1 then yields 61 = {I + l)/3. The limit x — > Xc in ()4.3|1 provides an 
expansion of the perimeter moments of the form 

where f3k,i = k(j) — 61. For the exponents describing the growth of the corresponding series 
coefficients a„ = [q'^]hk (q) in ()7.H) . (where [x"']5f(a;) denotes the coefficient of x" in the 
expansion of the function g{x)), it follows that a = k(j) and ctj = 6'j. We remark that the 
number of coefficients Mq used in the fit is much smaller than that for the area moments 
|15j . where 8 < Mq < 12. Apparently, the convergence of the perimeter moments to the 
asymptotic regime is quite slow. We were initially concerned that this significantly slower 
convergence was indicative of some feature of the scaling function we had overlooked. We 
were reassured that that is not the case, by performing the same analysis mutatis mutandis 
of the perimeter moment amplitudes for the (exactly solvable) model of staircase polygons. 
Precisely the same phenomenon was observed there, and in that case the scaling form has 
been proved [19J. 

For given M, the amplitude estimates {di} of (j7.2j) display cyclic fluctuations in A^. In 
order to enhance convergence, we considered only every r-th data value, i.e., we determined 
the coefficients di using sets of equations parametrised by n = N — r{M + l), N—rM, . . . ,N 
in fj7.2j) . where r = 2 for the square and hexagonal lattices, and r = 3 for the triangular 
lattice. The results of the fit are shown in Table |21 

For the coefficients Cq and ei, the scaling assumption leads to a prediction in terms of 
/o and /i from ()3.5|) . We get from ()3.5|) . on the square lattice, cq = —0.3941877(3) and 
ei = —2.575656(2). This agrees, within numerical accuracy, with the estimates obtained in 
Table 121 Similarly, for the hexagonal lattice, the estimates are consistent with the scaling 
function predictions Cq = —0.679256(2) and ei = —5.35661(1). For the triangular lattice, 
we get Co = —0.476162(2) and ei = —2.95663(2), which is again consistent with the result 
in Tabled 

It is often useful to check the behaviour of the amplitude estimates by plotting the 
results for the leading amplitude vs. 1/n. In Fig. ^ we have done so for the amplitude eo 
for the square, hexagonal and triangular lattices (the straight lines are the estimates given 
above). In the left panels we plot the estimates obtained with M ranging from 1 to 4 while 



Amplitude 


Exact value 


Square 


Hexagonal 


Triangular 


eo 


unknown 


-0.3942(2) 


-0.6790(4) 


-0.476(1) 


ei 


unknown 


-2.576(1) 


-5.356(2) 


-2.95(1) 


626^^60 


-0.29052826 


-0.29052(2) 


-0.29054(5) 


-0.2906(1) 


eaer^e^ 


0.1396478 


0.13967(4) 


0.13962(4) 


0.13965(3) 


—4 ■? 

6461 eg 


-0.0754402 


-0.07545(9) 


-0.07542(8) 


-0.0754(1) 


eser'e^ 


0.042564 


0.04259(7) 


0.0426(1) 


0.0426(1) 


e,e-,'el 


-0.02448 


-0.02451(7) 


-0.02446(8) 


-0.0245(2) 


6761 ^e^ 


0.0142143 


0.01423(8) 


0.01423(9) 


0.0144(3) 


ege^^e^ 


-0.008292 


-0.00829(9) 


-0.00831(6) 


-0.0082(4) 


eger^e^ 


0.0048499 


0.0048(1) 


0.00486(8) 


0.0048(4) 


eioer^°e^ 


-0.0028406 


-0.0028(2) 


-0.00284(7) 


-0.0029(3) 



Table 2: Bicritical perimeter moment amplitudes of rooted self-avoiding polygons 



the right panels give a closer look at the best converged sequences of amplitude estimates. 
In each case we use fits with a = —1/3 and ai = [i + l)/3 and as discussed above we have 
tried to minimise cyclic fluctuations. We observe that fits with M = 1 display pronounced 
curvature indicating that using just 1 sub-leading term gives an insufficient approximation. 
For the square and triangular cases the fits with M = 4 are marred by large fluctuations 
and are not very useful. The remaining fits clearly yield estimates for eo fully consistent 
with the precise values obtained above using the estimates for /o and /i. We notice that 
as more terms are added to the fits the estimates exhibits less curvature and that the slope 
become less steep (this is particularly so in the hexagonal case). This is evidence that we 
are indeed fitting to the correct asymptotic form. 

We finally considered the amplitude combinations e^e^ Cg"^, where 2 < A; < 10. They 
were estimated from the ratios 
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00 



(4.6) 



We extracted the amplitudes by a direct fit to the expected asymptotic form, as explained 
in the appendix. As argued above, we fitted with exponents of the form a, = {i + l)/3 for 
1 < i < Mq, where 2 < Mq < 4. The result is shown in Table El The prediction of the 
amplitude combinations appears to be correct, within numerical accuracy. 
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Figure 1: Estimates of the amplitude eo vs. 1/n for the square lattice (top panels) 
hexagonal lattice (middle panels) and triangular lattice (bottom panels). 



5 Unrooted self-avoiding polygons 

The kth perimeter moment of rooted self-avoiding polygons h^ (q) is related to the {k+l)th 
perimeter moment hk+i{q) of unrooted self-avoiding polygons. We have, for A; > 0, 



^i-'fe'-^^«'"(-" 



x—-G[x, q) 



k\ dx'^ V dx 



k hk{q) -{k + l) Xc hk+i{q) ~ -{k + 1) Xc hk+i{q) {q -^ I] 



(5.1) 



The last relation follows with the exponents l3k = 2k /3 — 1/3 of hj^ (q) given in ()3.2|) . It 
follows from ()5.1|) that, for A; > 0, the singular behaviour of hk{q) is determined by the 
singular behaviour of hf!^_i{q). So we have all the moments of unrooted SAP except the 
zeroth moment. It thus remains to analyse the moment ho{q). 

Extrapolating the values of /^^ to A; = —1 gives /3_i = —1, so we expect a singularity of 
the form 

hl--'\q) ~ Ail - q) log(l -q) (q ^ 1), (5.2) 

with some amplitude A > 0. This behaviour was tested by a direct fit to the expected 
asymptotic form. Using the notation of the appendix, we expect an exponent a = —1 and 
choose ai = i for 1 < i < Mq, where 2 < Mq < 4 for stable approximation schemes. 

For the square lattice, we find Cq = 0.026527(6). The numerical analysis yields Cq = 
0.026527(3) on the hexagonal lattice, and Cq = 0.05306(5) on the triangular lattice. This 
suggests a universal law of the form 

y'pm.na;™^^ 2 (n^oo), (5.3) 

m 

where a = 2 for the square lattice and the hexagonal lattice, and a = 1 for the triangular 
lattice. We have l/(127r) = 0.0265258... and l/(67r) = 0.0530516..., which is, within error 
bars, in agreement with the estimates obtained above. 

If one accepts the predicted scaling form (jl.ip and scaling function (j2.4|) of rooted self- 
avoiding polygons, (and we believe that we have provided compelling numerical evidence 
to do so), then the scaling function of (unrooted) self-avoiding polygons is determined by 
integration, 

G(-"^)(x, g) ~ (1 - q)F (j^^—^^ + C{q), (x, q) ^ (x„ 1), (5.4) 



'l-qY/\ 
where C{q) is a "constant" of integration, given below, and F{s) is given by 



'i^H^W) 



F(.) = -^logAi Uf . , (5.5) 
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where the coefficients /o and /i are given in Table ^ and Xc is given in the introduction. 
Alternatively, F{s) can be expressed in terms of the amplitude Eq given in Tabled by 

F{s) = -— log Ai (—{2EoY/h^ . (5.6) 

na \Xc J 

The term C{q) incorporates the singular behaviour of the bicritical area generating function 
ho{q) = G{xc, q) as q approaches unity, since the limit x — ^ Xc from the other term yields 
a vanishing contribution. The function C{q) is thus given by 

C{q) = ^{l-q)\og{l-q). (5.7) 

The scaling form (j5.4|) . (j5.6p and (j5.7j) refines the prediction given previously in [22^ by the 
addition of the term (j5.7|) . 

6 Crossover to branched polymer phase 

Let Xc{q) denote the radius of convergence of the perimeter and area generating function 
G^^\x,q) of rooted SAP for g < 1 fixed, as g —* 1. The scaling function prediction (jl.ip 
leads to a prediction of the slope of the critical line Xc{q), see also J21I2S|- The slope of the 
critical line is determined by the ffist singularity Sc of the scaling function H^^' (s) on the 
negative real axis. More precisely, for g < 1 fixed, the argument s = (xc — x)/(l — g)"^ of the 
scaling function is negative for x > Xc, attaining its singular value Sc for x = Xc{q) > Xc- 
We thus expect asymptotically 

x,(g)~x,-s,(l-g)* (q^l). (6.1) 

For the particular scaling function ()2.4|) . the point Sc is given by 

2 

/o 



, , , „, -2.3381074104... (6.2) 

From the values of Tabled we obtain s^ = -0.2608637(5), -0.2161405(20), and s^ = 
—0.274509(2) for the square, hexagonal, and triangular lattices respectively. Note that, 
due to the particular form of the scaling function, the same behaviour applies to unrooted 
self-avoiding polygons. 

Figure 121 displays a log-log plot of Xc{q) — Xc versus 1 — g. In each plot a straight line 
corresponding to the expected form 

Xc(g)-a;, = -s,(l-g)2/3 (6.3) 

is given. We get reasonable agreement with the predicted form. The estimates were 
obtained using third order differential approximants ^I], with the degree of inhomogeneous 
polynomial ranging from 5 to 15, and the requirement that averages must include at least 
85% of the approximants. For this part of our study we calculated the numbers Pm,n for 
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Figure 2: Plots of Xc{q) — Xc versus 1 — q for the square lattice (top), hexagonal lattice 
(middle) and triangular lattice (bottom). The solid line has the predicted slope 2/3. 
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m < 90 (square lattice), m < 134 (hexagonal lattice) and m < 50 (triangular lattice), and 
all relevant sizes of the area. 

If the scaling function H^'^\s) ~ a{s — Sc)^ about s = Sc, it follows that the singular 
part of G^''\x, q) behaves as G'('")'(**"'3)(a;, q) ~ a(l — q)~^'^{xc{q) — XcY for x -^ Xc{q). The 
scaling function ()2.4|1 has a simple pole at Sc so 

Qir),{s^ng)^^^ g) ~ a(l - g)^/^(x,(g) - x,)-\ (6.4) 

The differential approximant analysis typically yields reasonably accurate exponent esti- 
mates for q < 0.995 and confirms that G^^\x, q) has a simple pole as x — ;> Xc{q)- For exam- 
ple the analysis of the square lattice series yields g = —1.0(3) at g = 0.995, g = —1.03(2) 
at g = 0.99 and g = —0.999(2) at g = 0.95. For g closer to 1 the exponent estimates 
became unreliable in that the errors bars were as large as the estimates, e.g. at g = 0.997 
we found g = —1.6 ± 1.7. We also tried to calculate the amplitude a(g) ~ a(l — g)^^^, but 
unfortunately we could not get accurate estimates for g > 0.9, and so have been unable to 
numerically confirm the predicted behaviour. 

Finally, we checked the behaviour of the phase boundary Xc{q) as g — > 0. The coefficient 
of x"^ in G^^'\x,q) is a polynomial in g and as g — »• it becomes completely dominated 
by the term of lowest degree in g. We thus have to examine the behaviour of the terms 
flm = Pm.rimin' whcrc Pm,n^in ^^ ^^^ numbcr of polygons with perimeter m having the min- 
imal possible area n^nin- Clearly, the polygon formed by making a linear chain of unit 
cells contributes to a^- Unit cells on the square, hexagonal and triangular lattices have 
perimeter 4, 6 and 3, respectively, and it thus follows that rijam — ''^/2 (square), m/4 
(hexagonal) and m (triangular). On the square lattice it has been shown 9\ that a^ grows 
exponentially with m. am is bounded from above by the number of site trees of size n^i^ on 
the dual lattice and from below by the number of minimally spanning polyominoes of size 
nmin (a minimally spanning polyomino is a polyomino spanning a h x w rectangle having 
size h + w + 1). Similar arguments apply to the other lattices. So we form the gener- 
ating function S{t) = ^^ amt"^ and using differential approximants find that S{t) has a 
singularity at tc = 0.5189688(2) on the square lattice, tc = 0.6986253(5) on the hexagonal 
lattice, and tc = 0.346530(1) on the triangular lattice. Since t"^ = x™'g"'™'" it follows that 
^c{q) ~ tc/q^, with b = 1/2, 1/4 and 1 for the square, hexagonal and triangular lattice, 
respectively. Figure El shows a log-log plot of Xc{q) versus g. In each plot a straight line 
corresponding to the expected form Xc{q) = tc/q'' is also shown. 

7 Conclusions 

We have analysed bicritical perimeter moments of self-avoiding polygons using data ob- 
tained from exact enumeration on the square, hexagonal and triangular lattices. This 
yields a new check of the earlier scaling function conjecture for self-avoiding polygons. The 
numerical analysis supports the crossover behaviour of the critical line to the branched 
polymer regime. Whereas we find the scaling function conjecture for rooted self-avoiding 
polygons satisfied, it can be valid for unrooted self-avoiding polygons only in modified form. 
By analysing the bicritical area generating function, we suggest a modification by an ad- 
ditional term with an apparently universal amplitude, see (j5.4p . (j5.6p and (j5.7j) . It would 
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Figure 3: Plots of Xc{q) versus q for the square lattice (top), hexagonal lattice (middle) 
and triangular lattice (bottom). The solid line has the predicted form Xc{q) = tc/q^ with 
slope b = 1/2, 1/4 and 1, respectively. 

14 



be interesting to consider whether its value can be justified by field theoretical arguments; 
compare the related investigation [3]. 
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Appendix: Numerical methods 

We numerically analyse sequences a„ by a direct fit to the expected asymptotic form. 
Similar applications of this method can be found in |141 [T3] . The sequence a„ is assumed 
to behave asymptotically as 

n„a-l /" , Ci , C2 , , Cm 



a„~u"n""McoH \ h...H in -^oo), (7.1) 

with constants yU > 0, Cj 7^ for z = 0, 1, . . . , M, and exponents a and ai for z = 1, . . . , M, 
where ctj+i > a^ for i = l,...,M — 1. Estimates of the constant /i and the exponent a 
can be obtained by, e.g., the method of differential approximants pi]. We have /x = 1 in 
our examples. Often the sequence (aj)fl^ is unknown, but there are predictions for the 
numbers ctj. The validity of a prediction can be tested employing the following procedure: 
We perform a direct fit to the expected asymptotic form, i.e., we solve the linear system 



a„ = /i"n"-Mcio + ^ + ^ + --- + ^ (n = Ar-(M + l),...,Ar) (7.2) 

for the M + 1 unknowns di. If the assumption of the asymptotic form ()7.H1 is correct, then 
the numbers di = di{N, M) will satisfy 

di{N,M)^Ci [N-^oo). (7.3) 

Generally, if the wrong sequence has been chosen, the sequence of coefficients di{N,M), 
for increasing values of N and fixed M, diverges either to infinity or converges to zero; but 
if the correct sequence has been chosen, convergence is usually rapid and obvious. 

Let us fix i in the following. Estimates of the amplitudes Cj are obtained in the following 
way. For A^ large enough, the numbers di{N, M) display approximately linear variation in 
l/AT, 

d,{N, M) ~ q(M) + ^^ (AT ^ 00). (7.4) 
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The numbers Ci{M) are obtained by a linear least squares fit of di{N, M) for large values 
of A^. The larger M, the larger A^ has to be taken in order to reach the asymptotic regime 
()7.4|1 . Since, however, N < Nq for given series data, di{N,M) is close to the asymptotic 
regime only for values of M < Mq, where Mq has to be extracted from series analysis. We 
choose 1 < mo < Mq such that |ri(mo)| is minimal. We estimate q by Ci(mo) and estimate 
the error by the spread among different values Cj(m) for 1 < m < Mq. 

The amplitudes Ci in ()7.1|) are related to the critical amplitudes of the corresponding 
generating functions. If a; ^ — Nq, the singular amplitude B of the corresponding generating 
function 

oo 

A{x) = ^a„x" ~ A(-^)(x) + _ {fix ^ 1) (7.5) 

n=0 ^^ ^' 

is related to Cq via B = coT{a)fi~". 
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